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Abstract 

The physics of disordered alloys, such as typically the well known case of 
CeNil-xCux alloys, showing an interplay among the Kondo effect, the spin glass 
state and a magnetic order, has been studied firstly within an average description 
like in the Sherrington-Kirkpatrick model. Recently, a theoretical model [1] in- 
volving a more local description of the intersite interaction has been proposed to 
describe the phase diagram of CeNil-xCux. This alloy is an example of the com- 
plex interplay between Kondo effect and frustration in which there is in particular 
the onset of a cluster-glass state. Although the model given in Ref. 1 1 1 has repro- 
duced the different phases relatively well, it is not able to describe the cluster-glass 
state. We study here the competition between the Kondo effect and a cluster glass 
phase within a Kondo Lattice model with an inter-cluster random Gaussian inter- 
action. The inter-cluster term is treated within the cluster mean-field theory for 
spin glasses |2|, while, inside the cluster, an exact diagonalisation is performed 
including inter-site ferromagnetic and intra-site Kondo interactions. The cluster 
glass order parameters and the Kondo correlation function are obtained for dif- 
ferent values of the cluster size, the intra-cluster ferromagnetic coupling and the 
Kondo intra-site coupling. We obtain, for instance, that the increase of the Kondo 
coupling tends to destroy the cluster glass phase. 

1 Introduction 

The properties of many cerium or uranium compounds are well described by the Kondo- 
lattice model, in which there is a strong competition between the Kondo effect on each 
site and the Ruderman-Kittel-Yosida-Kasuya (RKKY) interaction between magnetic 
atoms at different sites. On the other hand, it is now quite clear that the interplay be- 
tween disorder and electronic correlations in these systems produces a new physics 
such as the presence of Non-Fermi liquid behaviour O SI [5) or percolative process in 
the magnetic states similar to manganites [6|. 

An example of the percolative scenario can be found in CeNi\- x Cu x Q. This 
particular alloy presents a quite complex phase diagram as long as the doping of Ni 
increases. In the range of doping 0.6 < x < 0.3, the /iSR results indicate formation of 
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clusters below a characteristic temperature T*. On the other hand, the ac-susceptibility 
(Xac) shows the presence of a glassy ordering below a freezing temperature 7V (Tf < 
T*). In this sense, this particular state can be caracterized as a Cluster Glass (CG) state. 
Finally, it has been found an inhomogeneous ferromagnetic (IFM) order from neutron 
difraction at much lower temperature. However, there is no clear indication of a Curie 
temperature T c from Xac an d C p measurements. Therefore, it is possible to speculate 
that the evolution from CG state to the IFM order can be obtained by the percolation 
of the frozen clusters. One important experimental evidence supporting such scenario 
has been the behaviour of the hysteresis loops for CeMo.6C«o.4 at T = 100 mK which 
display discrete jumps of the magnetization. 

From the theoretical point of view, a model has been proposed to explain the pres- 
ence of frustration in the global phase diagram of CeNi\- x Cu x . The model is a Kondo 
Lattice with an additional Ising intersite interaction between localized spins called here 
Kondo-Ising Lattice (KIL) model [9 1. The important point is that disorder can be intro- 
duced in the KIL model by choosing, for instance, the coupling between the local- 
ized spins as a Gaussian random variable (see for example, the Sherrington-Kirkpatrick 
(SK) model [ 10 1). The results have shown that it is possible to construct a mean field 
solution of the KIL model, where it is found a spin glass SG solution as well as a Kondo 
regime [9 |. This approach has been extended to include a ferromagnetic (FM) solution 
[ 1 1 1 by displacing the Gaussian distribuition from the origin to Jq. This procedure has 
allowed to introduce the usual magnetization as a new order parameter. As a result, a 
global phase diagram temperature versus Jk (the strength of Kondo coupling) has been 
obtained displaying a SG phase, an additional FM one and a Kondo regime. However, 
the sequence of phase transitions when the temperature is decreased (for a constant Jk) 
is not in agreement to the experimental findings. On the contrary, the FM solution ap- 
pears at higher temperature than the SG one. Moreover, there is no percolation process. 
The results also show a conventional phase transitions in which there is, for instance, 
a clear Curie temperature T c , It is also important to remark that the usual FM order 
parameter included in this particular approach is not able to capture the complexity of 
the experimental IFM ordering. 

Recently, a new approach has been proposed replacing the random Gaussian 
in Eq. (OQ) for random site model in which /y — £^ =1 §j where %f are random 
variables which follow the distribution P{£,f ) — 1 /28^n — 1 /25cm E2- ^ n f act ' 
this particular choice of /y can allow the interpolation from weak to strong frustration 
regimes. This model improves the previous SK-based model in two directions: first 
it gives a better possible description of the experimental IFM ordering and second it 
yelds a disordered ferromagnetic phase below the spin glass one in better agreement 
with the experiment. However, it is still necessary further theoretical improvements to 
describe the magnetic clusters in CeNi\- x Cu x . 

The presence of the cluster glass is a clear indicative that the frustration present in 
the intermediated doping of CeNi\- x Cu x can not be described by a conventional spin 
glass. One possible improvement to the original KIL model would be to reformulate the 
intersite random interaction using cluster of spins instead of canonical spins. In fact, in 
a earlier work [2], the classical cluster glass problem has been studied in a mean field 
level. The model used is composed, basically, by a intracluster ferromagnetic coupling 
and a intercluster Gaussian random coupling. This kind of approach seems adequate 
to be implemented to study the competition between Kondo effect and cluster glass 
within the approach of Ref. |9|. Therefore, cluster of spins could be introduced in the 
original KIL model 13 by replacing the intersite random Gaussian term by intracluster 
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and intercluster terms similar to Ref. J2)- 

The aim of the present work is to study competition between Kondo effect and 
cluster glass. In view of the discussion in the previous paragraph, we use the following 
model (called here Kondo Lattice Cluster Glass (KLCG)) to accomplish that study: 



N a n s N c i n, N d 

o=li=l <J=TI a=\i<j a i.j 



(1) 



a=\i=\ a<b 



where intercluster coupling J ab is a random variable as the Sherrington-Kirkpatrick 
model 

nW^e-^^^ (2) 

with N = N c i * n s , where N c [ and n s are the number of cluster and the number of spin 
in each cluster, respectively. The hopping tf" is only inside the cluster. The indices 

(a,b) refer to clusters while indicates spins inside a cluster. So, S^ a — f} a ^fia[, 

s 7a = 4j/faT and 

<=i<t=U '=i 

SI = E Oflaafica (4) 

The problem can be treated within the formalism of integral functional where the 
spin operators are given by bilinear combinations of Grassmann fields. It should be 
remarked that disorder is introduced only in the intercluster random Gaussian interac- 
tion J a f,. Nonetheless, the thermodynamics for this particular situation is also obtained 
using the replica method. So, 

1 (Z")j , - 1 
BF = lim . (5) 

N n^o n 

The procedure to deal with the disorder follows closely the the usual fermionic 
spin glass (see for instance lfl3l ). The problem is treated within the static approxi- 
mation with order parameters q a p (a ^ j3) and q aa being introduced by a Hubbard- 
Stratonovich transformation. Then, replica symmetry is assumed q = q a p and p = q aa . 
The main difference is that, in the present approach, these order parameters describe 
a glassy ordering among cluster instead of canonical spins. The details of the calcula- 
tions will be shown elsewhere lfT4l . 
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Figure 1: Spin glass order parameter and correlaction function (SjSj) versus h/J f° r 
n s = 2 and Jq/J — 1. The dashed lines are results at temperature T = 0.507 and the full 
lines correspond to results at T = 0.757. 



and 



2% ff = |D(<p»D(v/» e A ^ 
That is equivalent to diagonalize the following Hamiltonian: 



(7) 



(=i 



'■] 



(8) 



'=1 ' i<j .7=1 



where h(p,q) is given as 



h(p,q) = pjy/2(p-q)Z + JTqz. 



(9) 
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Then, in the effective problem, h(p,q) appears as a random external field applied in 
the cluster which depends on the clusters glass order parameters q and p. Therefore, to 
obtain any information in this problem, it is necessary to diagonalize the Hamiltonian 
H e ff given in Eq. ^ and, simultaneously, to solve the saddle point equations for q and 
P- 



5 
4 

cr 3 
2 
1 




10 



15 
J,,/J 



~~\ r 

T=0.75J - 
T=0.50J - 



20 



25 



30 



-0.25 
-0.2 

, A _ -0.15 

v 

-0.1 

-0.05 



in 




30 



Figure 2: Spin glass order parameter and correlaction function (SJs?) versus J^/J for 
n s = 3 and Jq/J = 1. The dashed lines are results at temperature T = 0.507 and the full 
lines correspond to results at T = 0.75/. 



The behaviour of the cluster glass order parameters q and the Kondo correlation func- 
tion < Sfs^ > are displayed in Figs. 1-2 as a function of the Kondo coupling Jk for two 
values of temperature T while the intracluster ferromagnetic coupling is kept constant 
Jo = / (/is defined in Eq. (f2|l). The size of the clusters also assumes two values, n s = 2 
and 3. In Fig. 1, the results for n s = 2 show that, when increases, the cluster glass 
order parameters q decreases. For T = 0.757 and T = 0.5/, the cluster glass phase 
are destroyed for Jg w 11/ and Jk « 13/, respectively. While q decreases, < Sjsj > 
enhances which means that the Kondo effect inside the cluster becomes increasingly 
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important. The combined behaviour of q and < iS?s? > would indicate that increase 
of the Kondo effect inside the cluster is related to the destruction of the cluster glass 
phase. For ns = 3, the scenario described previously is preserved with q decreasing and 
< Sp 7 : > increasing. However, q vanishes for larger values of Jk as compared with the 
case n s = 2, particularly, for T — 0.57. That result would suggest that the cluster glass 
phase becomes more robust with the increase of n s 

To conclude, in the present work we have studied the competition between cluster 
glass phase and Kondo effect using the Kondo Lattice Cluster Glass model defined 
in Eq. (Q~|). The intercluster disorder problem is treated using the usual mean field 
procedure for fermionic spin glasses. Therefore, the original problem is transformed in 
an effective one in which there is a random external field applied on the cluster. Finally, 
it is used exactly diagonalization to solve the cluster. The results indicate that, when 
Jk increases, the Kondo correlation function < S^s\ > also increases. Simultaneously, 
the cluster glass phase is destroyed. These results suggest that this approach could be 
used to study the behaviour of the CeNi\- x Cu x . However, as discussed previously, the 
frustration in cited physical system can not be described in terms of a random Gaussian 
coupling. It would be better described by a coupling used in Ref. [ 1 1. This approach is 
under current investigation. 
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